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ABSTRACT 

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and 
an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from 
Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of 
the AGN over two decades in energy (^0.5-100 keV). Previous X-ray observations suggested that 
the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. 
However, the lack of high-quality > 10 keV observations, together with the presence of a nearby 
X-ray luminous source, NGC 5643 X-I, had left significant uncertainties in the characterization of 
the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved 
above 10 keV for the first time and allows a direct measurement of the absorbing column density 
toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of A^h ^ 

5 X 10^® cm“^. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting 
models is T 2 -io,int = (0.8-1.7) x 10®^ erg s“®, consistent with that predicted from multiwavelength 
intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X~l, and show that 
it exhibits evidence for a spectral cut-off at energy, A ^ 10 keV, similar to that seen in other ULXs 
observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 
keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1. 

Subject headings: galaxies: active — galaxies: nuclei — techniques: spectroscopic — X-rays: galaxies 
— X-rays: individual (NGC 5643, NGC 5643-Xl) 


1. INTRODUCTION 

Compton-thick AGN (CTAGN) are expected to consti¬ 
tute a significant fraction of the overall AGN population 
in the local universe, accounting for ^20 -30% of AGN 
according to inultiwavele n gth studies (e.g . iRisaliti et al.l 
119991: iBurlon et al.ll201ll : Goulding et al'l [201111 . Many 
studies also predict that CTAGN provide a substantial 
contribution to the cosmic X-ray background (CXB), 
responsible for 10-25 % of the flux at the peak en- 
ergy, ^ 30 keV (e.g. iGilli et all 120071: iTreister et al.l 


2009t iDraoer fc Ballantvnd 120101 : lAkvlas et al.l 120121 : 
Ueda et al.ll2014D . Yet, their census is still far from com¬ 
plete. The high line-of-sight column density in CTAGN 
(-^H > 1/<7 t = 1-5 X 10^® cm“^, where ax is the 
Thomson scattering constant), generally attributed to 
the parsec-scale circumnuclear torus of AGN unification 
schemes, as well as larger-scale molecular clouds and 
dust lanes, results in severe attenuation of the direct X- 
ray emission from CTAGN at energies below 10 keV. 
This is why observations at higher energies are needed 
to probe this direct component and provide unambigu- 
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ous identification of CTAGN. However, even at A > 10 
keV, the remaining flux that we observe in the most 
extreme CTAGN will be from photons scattered or re¬ 
flected from the back-side of the torus, and will com- 
prise just a few perc e nt of the intrinsic AGN power (e.g. 
Iwasawa et al.l 119971 : IMatt et HI 120001: iBalokovic et alJ 
2014) . This makes the identification and characteriza¬ 
tion of CTAGN a challenging task. 

To date, only ^20 AGN within 200 Mpc have been 
confirmed as C T based upon detailed X-ray spectral 
characterization ([ Della C eca et al.ll2008l : iGoulding et al.l 
120121: iGandhi et al.l 1201^ This corresponds to a frac¬ 
tion of <C 1% of the total AGN population expected 
within that volume, suggesting that the vast majority of 
CTAGN are yet to be found even in the local universe. 
These hona-fide CTAGN were unambiguously identified 
based upon a detection at energies above 10 keV and the 
presence of a Fe Ka line at 6.4 keV with high equiva¬ 
lent width, EW > 1 keV. The identification and char¬ 
acterization of all the CTAGN is important in order to 
form an accurate census of accretion in the local universe, 
since much of the growth of supermassive black holes is 
thought to occur in such heavily obscu red phases (e.g. 
lFabianiri999l : [Alexander fc Hickoxll201^ . Therefore, an 
accurate local benchmark is important for extrapolating 
the results to higher redshifts. 

We have started a program to study a complete, 
volume-limited {D < 15 Mpc), mid-infrared (MIR) se - 
lected AGN sample from IGoulding fc AlexandSl (j2009H . 
with the main goal of constraining the population of 
CTAGN and the A^h distribution of AGN in the lo¬ 
cal universe. CTAGN candidates from the sample were 
identified using multiwavelength selections, such as X- 
ray spectroscopy, and intrinsic 2-10 keV luminosity in¬ 
dicators from high spatial resolution MIR 12/mi contin¬ 


uum and [Oiii]A 5007A line luminosity corrected for the 
Balmer decrement. One of the candidates that stands 
out in the sample as being CT based on these analy¬ 
ses is NGC 5643. NGC 5643 is a nearby face-on (z « 
30° (Q SAB (rs)C galaxy hosting a low-luminosity Seyfert 
2 nucleus (|PhilliDS et ahlfTo^ . It has a redshift oi z = 
0.0040, corresponding to a metric/proper distance of D 
= 13. 9 Mpc under the assum ption of the iMould et al.l 
(|2000D cosmic attractor model ((Sanders et al.ll2003D I4 
NGC 5643 features a compact radio core with two- 
sided . kiloparsec-scale lobes in an east-west orienta¬ 
tion ((Morris et al.lll98,^ . and a cospatial one-sided Ha 
and [Oiii] emission line region e xtending east-ward o f 
the nucleus for at least 1.8 kpc ((Simpson et aLlllQ^ . 
Despite the intense star formation episodes occuring 
in the spiral arm, MIR diagnostics suggest that the 
AGN still d ominates the overa ll IR (8-1000p,m) en¬ 
ergy budget ((Genzel et al.lll99^ . Comparisons of op¬ 
tical spectra with synthesis models, however, are con- 
sistent with a “starburst/Seyfe rt 2 composite” spec- 
trum (iC id Fern andes et ^120011) . Using Bry emission, 

I Davies et al.l ( 20141 ) found no on-going star formation ac- 


^ The host galaxy inclination was obtained from the HyperLeda 
website (http://leda.univ-lyonl.fr/). 

^ IMould et al.l II2000I 1 adjusts heliocentric redshifts to the cen¬ 
troid of the Local Group, taking into account the gravitational 
attraction towards the Virgo Cluster, the Great Attractor and the 
Shapley supercluster. 


tivity in the nucleus, although the possibility of a recent 
(terminated) starburst cannot be excluded. This source 
also shows spatially resolved molecular ga s flowing out 
from the AGN at a rate of 10 Mq yr“^ ((Davies et al.l 
l2014|) . Although water maser emission associated with 
the AGN has been detected in NGC 5643, a correspond¬ 
ing spatially resolved map, which would allow for a di¬ 
rect measurement of the s upermassive black ho le mass 
(Mbh), is not yet available (iGreenhill et al.l[2003D . How¬ 
ever, an indirect Mbh measurement from the galaxy stel¬ 
lar velocity dispersion (cr*) provides an estimated black 


hole mass of Mbh = 10®-^ Mq ((Goulding et al.l[2010D Pl 

In X-rays, NGC 564 3 has been obse rved by ASCA 
(iGuainazzi et al.l 120041). BeppoSAX (iMaiolino et al.1 
1 1 9981 ). ROSAT (Guainazzi et HI I2004D . Chan¬ 
dra (iB ianchi et al. I2006D and XMM-Newton 
((Guainazzi et al.l 12004 lMa,tt et al.l I2013D . Dra¬ 

matic spectral changes were observed between the 
XMM-Newton observation carried out in 2003, and the 
ASCA and BeppoSAX observations performed earlier. 
However, the point spread function (PSF) of ASCA 
and BeppoSAX were not sufficient to separate the 
emission of the nucleus from that of a nearby X-ray 
source (at an angular separation of ^50"), NGC 5643 
X-1, which was found to be very bright at the time 
of the XMM-Newton observation (G04). Therefore, 
it remained unclear which source was responsible for 
the spectral variability observed. Comparisons of the 
XMM-Newton observations in 2003 and 2009 showed 
that there is no significant v ariation in the spectrum of 
the AGN ((Matt et al.l 1201,K : hereafter M13). However, 
the off-nuclear source was found to be more than a 
factor of two fainter in flux in 2009 (M13) than in 2003 
(G04). 

The Chandra image of the AGN shows that the soft X- 
ray emission {E <2 keV) of th e nucleus is spatially corre¬ 
lated with the [Oiii] emission ((Bianchi et al.l[2(if)ffl . con¬ 
sistent with what is commonly observed for many Seyfert 
2 galaxies. The dominant power source of this soft X-ray 
emission appears to be photoionization from the AGN, 
although it is still unclear how much collisional ionization 
contributes to the overall X-ray emission (M13). Above 
2 keV, the X-ray spectra from these various observations 
show indications of the nucleus being absorbed by CT 
material. The evidence for this are the detection of a 
prominent Fe Ka line (EW > 1 keV) and a flat photon in¬ 
dex below 10 keV (r 2 _io < 1), which are characteristics 
of a reflection-dominated spectrum. Analysis of the low 
signal-to-noise ratio (S/N) BeppoSAX spectrum where 
NGC 5643 was detected only up to 10 keV, combined 
with an upper limit for the 15-100 keV band, suggested 
a tentative lower limit to the co lumn density of Nn > 
10^® cm“^ ((Maiolino et al.l[i998l) . 

The true nature of the off-nuclear X-ray source, NGC 
5643 X-1, is still uncertain. It is located in the outskirts 
of the host galaxy optical emission and identification of 
counterparts at other wavelengths has been ambiguous 
(G04). It is highly likely that the source is located inside 
the galaxy and therefore would be a powerful ultralumi- 


® The same Mbh estimate is obtained using an updated 
value of i7„ de ri ved fr om the [O iii]A 5007A emission line width 
from IG u et al.| II2006I) and the latest Mbh-c* correlation by 
IMcConnell Kr. Ma Il2ni3l ~). 



















































































NuSTAR Observations of NGC 5643. 


3 


nous X-ray source (ULX) with ^^ 2-10 ~ 1-7 x 10“^° erg 
s“^ based on the flux observed in 2003 (G04). This is 
comparable to the observed luminosity of the AGN itself; 
i.e., ^2-10 ~ 1.9 X 10^° erg s“^ (G04). 

In this paper, we present new NuSTAR observations of 
NGC 5643 in which the AGN and the off-nuclear source, 
NGC 5643 X-1, are clearly resolved and detected at hard 
X-ray energies (> 10 keV) for the first time. This allows 
us to provide the most accurate spectral analysis of the 
AGN to date. The aim of this paper is to characterize 
the broadband spectrum of the AGN by combining our 
NuSTAR data with existing data from Chandra, XMM- 
Newton and Swift-BAT. We also present the NuSTAR 
data for NGC 5643 X-1, which is detected above 10 keV 
for the first time. 

The paper is organized as follows: we describe details 
of the X-ray observations and data reduction of the AGN 
in Section 2, followed by the spectral fitting procedures 
and results in Section 3. In Section 4, we present the data 
analysis and results on NGC 5643 X-1. This is followed 
by a discussion in Section 5. The paper concludes with 
a summary in Section 6. 

2. OBSERVATIONS 

In this section, we describe the NuSTAR observations 
and data analysis procedures for the AGN. We also detail 
the archival Chandra, XMM-Newton and Swift data that 
were used to facilitate our broadband spectral analysis. 
Details of these observations are provided in Table 1, and 
the data reduction is described below. 

2.1. NuSTAR 

NuSTAR (Nuclea r Spectroscopic Telescope Array, 
iHarrison et al.l[2013f) . operating at 3-79 keV, is the first 
orbiting observatory with the ability to focus X-ray pho¬ 
tons at A > 10 keV. It provides a two orders of mag¬ 
nitude improvement in sensitivity and over an order of 
magnitude improvement in angular resolution with re¬ 
spect to previous hard X-ray orbiting observatories. It 
consists of two co-aligned focal plane modules, FPMA 
and FPMB, each covering the same 12' x 12' portion 
of the sky, and each comprised of four detectors placed 
in a 2 X 2 array. NuSTAR has an angular resolution 
of 18" FWHM with a half power diameter of ^1', and 
energy resolutions (FWHM) of 0.4 and 0.9 keV at 6 
and 60 keV, respectively. These capabilities make NuS¬ 
TAR an ideal instrument to characterize the spectral 
shape of heavily obscured AGN in the local universe (e.g. 
Balokovic et ai.l[20l4 iGandhi et alll20l4 IPuccetti et al.l 

20141 iBauer et al.ll2ni4lb 

NGG 5643 was observed twice by NuSTAR in 2014 
- an initial observation with a nominal exposure time 
of 22.5 ks taken in May, followed by an additional 19.7 
ks observation in June. The second observation was con¬ 
ducted to improve the S/N of NGG 5643 as it fell near the 
detector gap in the first observation. The data were pro¬ 
cessed with the NuSTAR Data Analysis Software (nus- 
TARDAS) vl.4.1 within HEASOFT v6.15.1 with GALDB 
v20150316. Galibrated and cleaned event files were pro¬ 
duced using the NUPIPELINE vO.4.3 script with standard 
filter flags. Spectra and response files were extracted us¬ 
ing the NUPRODUCTS vO.2.5 task. 

For each observation, the AGN spectra were extracted 
using a circular aperture region of 30"-radius centered 


on its peak emission. The background photon were col¬ 
lected from polygon-shaped regions lying on the same 
detector as the source for the June observation, and from 
two adjacent detectors for the May observation (since the 
source fell close to the detector gap in this observation). 
We excluded background photons that lay within a cir¬ 
cular region of ~ 70"-radius around the source to exclude 
emission from NGC 5643 X-1. Significant counts are de¬ 
tected up to 50 keV for the AGN. We show the combined 
FPMA-I-B images of the AGN and NGC 5643 X-1 in the 
3-8, 8-24 and 24-50 keV bands in Figure 1. 

2.2. XMM-Newton 

NGC 5643 was also observed on two occasions by 
XMM-Newton, for 10 ks in 2003 and 55 ks in 2009. In 
both cases, the observations were performed with the 
EPIC CCD cameras (PN, MOSl, MOS2), operated in 
full frame mode with the medium filter. The observations 
were discussed in detail in G04 and M13, respectively. 
For the AGN, we only extracted and used the spectra 
from the longer 2009 observation as it has the highest 
S/N. Data were reduced within SAS vl.2, screened for 
flaring particle background as described in M13. 

Source spectra were extracted within a 25"-radius 
aperture centered on the AGN for all cameras. Back¬ 
ground photons were selected from source-free circular 
regions of 100"-radius on the same chip as the source. 
Patterns 0 to 4 and 0 to 12 were used for the PN and 
MOS spectra, respectively. The spectra from MOSl and 
MOS2 were co-added as the data were consistent with 
each other (see Section 3 for further details). 

2.3. Chandra 

NGG 5643 has only been observed once by Chandra. 
This observation was conducted in 2004 with the AGIS-S 
detector with an exposure time of ~ 8 k s. The results 
of the observation were first published in IBianchi et al.l 
(1200611 . We reprocessed the data to create event files with 
updated calibration modifications using the GIAO v4.6 
pipeline following standard procedures. 

Source counts were extracted using the SPECEXTRACT 
task in GIAO from a circular region of 25"-radius cen¬ 
tered on the AGN to match the XMM-Newton extraction 
region. The background was extracted from a source-free 
30"-radius circular region close to the source. 

2.4. Swift 

The first NuSTAR observation in May 2014 was ac- 
companied by a ^ 5 ks Swift X-ray Telescope (XRT; 
iBurrows et al.l 1^00511 observation, which started approx¬ 
imately an hour after the NuSTAR observation began. 
The purpose of this observation was to provide simul¬ 
taneous coverage for the soft X-ray end of the spec¬ 
trum (< 3 keV) where the NuSTAR sensitivity drops off. 
The data were reduced using the xrtpipeline vO.13.0, 
which is part of the XRT Data Analysis Software (xRT- 
DAS) within HEASOFT. However, with only ~ 20 counts 
between 0.5-10 keV, the exposure is not long enough 
to provide additional constraints beyond those already 
obtained with NuSTAR, XMM-Newton and Chandra. 
Therefore, we only used these data to check for con¬ 
sistency with the XMM-Newton and Chandra data at 
0.5-10 keV (see Section 3). 
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Figure 1. NuSTAR and XMM-Newton images of the AGN (NGC 5643) and the ULX candidate (NGC 5643 X-1). The sources are 
circled in green with 30'' and 20"-radius apertures, respectively. Images are smoothed with a Gaussian function of radius 5 and 3 pixels, 
corresponding to 12.3" and 3.3" for NuSTAR and XMM-Newton, respectively. North is up and East is to the left in all images. Top: 
NuSTAR combined FPMA+B images in the 3-8, 8-24 and 24-50 keV band for the May 2014 observation (first row) and June 2014 
observation (second row). Blue dashed lines correspond to the gap between the NuSTAR detectors for the May observation. Bottom: 
XMM-Newton EPIC MOS 2 image (July 2009) in the 0.5-10 keV band shown to provide a clearer view of the two sources at lower energies. 


Table 1 

X-ray Observations log of NGC 5643 


instrument 

(1) 

OTilD 

(2) 

Date 

(3) 

Energy Band 
(keV) 

(4) 

Net Exposure time 
(ks) 

(5) 

Net count Hate 
(10-2 cts s-l) 
(6) 

Chandra ACIS-S 

5636 

2004-12-26 

0.5-8 

7:63 

1:96 

XMM-Newton PN/MOSl+2 

0601420101 

2009-07-25 

0.5-10 

45.4/53.4 

14.6/3.76 

Swjft-BAT 

- 

2004-2010 

14-100 

7340 

0.00223 

Swift-XRT 

00080731001 

2014-05-24 

0.5-10 

3.96 

0.472 

NuSTAR FPMA/FPMB 

60061362002 

2014-05-24 

3-50 

22.5/22.4 

2.44/1.94 

NuSTAR FPMA/FPMB 

60061362004 

2014-06-30 

3-50 

19.7/19.7 

2.32/2.36 


Note. — (1): List of observatories and instruments; (2): Observation identification number (obsID); (3): Observation UT start date; 
(4) Energy band in keV; (5) The net exposure time in ks; (6): Net count rate for the AGN in the given energy band in units of 10“^cts 
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The Burst Alert Telescope (BAT; iBarthelmv et al.l 
1200511 onboard Swift has been continuously monitoring 
the sky at 14-195 keV, producing images of a large num¬ 
ber of hard X-ray sources thanks to its wide field of 
view and large sky coverage. We used the stacked 70- 
month spectrum and its associated response file down¬ 
loaded from the Swift-BAT 70-month Har d X-ray Survey 
Source Catalog ((Baumgartner et ^1201311 to provide X- 
ray constraints above the NuSTAR bandQ NGC 5643 
is detected at the 5.4cr significance level with significant 
counts up to ~ 100 keV. 

3. BROADBAND SPECTRAL MODELING OF THE AGN 

We describe the broadband X-ray spectral analysis of 
the AGN in this section. The analysis was carried out 
using XSPEC vl2.8. 20 All spectra were binned to a mini¬ 
mum of 25 counts per bin to allow the use of statistics. 
Absorption through a fixed Galactic col umn along the 
line-of-sight, = 8.01 x 10^° cm“^ (iKalberla et al.l 
[200^, was included in all spectral fits using the xspec 
model “PHABS” , and solar abundances were assumed for 
all models. All errors are quoted at 90% conhdence, un¬ 
less stated otherwise. Details of the main results are sum¬ 
marized in Table 2, and the best fit spectra are shown in 
Figure 2. 

We began our spectral modelling with an examination 
of the NuSTAR data alone. We started by modeling the 
spectra of each of the NuSTAR observations in the 3- 
50 keV band using a simple absorbed power-law model 
with the column density fixed to the Galactic value. A 
prominent excess of counts just above 6 keV was ob¬ 
served in both spectra, suggesting the presence of Fe Ka 
emission. The parameters returned by the two spectra 
were consistent with each other, indicating that there are 
no significant differences between the two observations. 
Therefore, we co-added the spectra for each FPM using 
the ADDASCASPEC Script (the same test and procedures 
were also done for the two XMM-Newton MOS spectra). 

We then modeled the combined NuSTAR spectra in 
the 3-50 keV band with an absorbed power-law model 
(column density fixed to N^^^) and a Gaussian line at 
E ss 6.4 keV. The model measured a photon index and 
Fe Ka EW of F = 0.55 ± 0.07 and EWpe = 2.22^°;^^ 
keV, respectively, with a fit statistic of = US for 
76 degrees of freedom (d.o.f). The line energy is cen- 
tered at E = 6.36 dz 0.04 (statistical) ± 0.04 (systematic; 
iMadsen et al.l[20T5ll keV, consistent with neutral Fe Ka 
emission. The very flat F and large EWpe, which were 
also found in previous observations, are characteristic sig¬ 
natures of Compton-thick absorption and reflection from 
optically thick cold gas. 

We checked for variability in the flux of the AGN in 
the 0.5-10 keV band by comparing the XMM-Newton 
and Chandra data with the S'wift-XRT data. The fluxes 
of all spectra are consistent with each other, f 0 . 5 - 10.0 
10“^^ erg s“^ cm“^, suggesting that there are no sig- 
nihcant differences between the three observations. We 
also checked for variability at hard X-ray energies by 
constructing power spectra from the 3-50 keV NuSTAR 

^ The Swift-BAT 70-month Hard X-ray Survey Source Catalog 
is available online at http://swift.gsfc.nasa.gov/results/bs70mon/ 

® The XSPEC manual can be downloaded from 
http: //heasarc.gsfc.nasa.gov/xanadu/xspec/XspecManual.pdf 


lightcurves from the two observations. These were then 
compared to the power spectra expected for pure Pois¬ 
son noise and for the expected transmitted component 
variability observed in an unobscured AG N with a simi¬ 
lar bl ack hole mass and accretion rate (see lArevalo et al.l 
120141 for further details). No evidence for variability was 
found, however, given the low count rates, this test is not 
very sensitive. Given that the NuSTAR spectra from the 
two observations are consistent with each other, we as¬ 
sumed that the AGN has not varied. 

We therefore proceeded to fit the XMM-Newton and 
Chandra spectra combined with the NuSTAR and Swif- 
BAT spectra using more detailed physical models to bet¬ 
ter characterize the broadband spectrum of the AGN. 
We describe the details and resu lts of each model used: 
Mode l P (PEXRAV model bv iMagd ziarz fc ZdziarsO 
I1995H . T (torus model bv iBrigh tman fc Nandra 2011fl 
and M (MYTORUS model bv iMurnhv fc YaaoobI 200^ . 
in Section 3.1, 3.2 and 3.3, respectively. In addition to 
these models, we added extra components required to 
provide a good fit to the data, briefly described below. 

The soft energy (< 2 keV) part of the spectra covered 
by XMM-Newton and Chandra is dominated by emis¬ 
sion from photoionized material (M13). We modeled this 
emission based on M13 using a soft power-lav0 and 10 
Gaussian components (9 for Model M; see Section 3.3) to 
model the emission lines (refer to Table 1 in M13). We 
also added four and one more Gaussian component (s) to 
Model P and M, respectively, to model the fluorescence 
lines emitted at E > 2 keV which are not included in 
these models (see Section 3.1 and 3.3 respectively). 

Cross-calibration uncertainties between each obser¬ 
vatory with respect to NuSTAR were included as free 
parameters (constant, C). Initially, we left both the 
XMM-Newton EPIC-MOS and PN cameras constants 
free to vary. However, we found that their values are 
consistent with each other and decided to tie them 
together in the final fittings of all models. 

Our three models can be described as follows: 

Model P = CONSTANT X PHABS X (POW -|- 

PEXRAV -I- 14 X ZGAUSS), (1) 

Model T = CONSTANT X PHABS X (POW -|- 

TORUS -I- 10 X ZGAUSS). (2) 

Model M = CONSTANT X PHABS X (POW -|- 

ZPOW X MYTZ -I- MYTS -|- MYTL -|- 

10 X ZGAUSS), (3) 


® The soft power-law used to model the unresolved < 2 keV 
emission also includes the scattered power-law component which 
is commonly used to simulate the AGN emission scattered into 
our-line-of sight by diffuse hot gas. 
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Figure 2. Best-fit models to the combined NuSTAR, Chandra, XMM-Newton and Swift-BAT data of the AGN - Model P (top), Model 
T (bottom left), and Model M (bottom right). Model P and T were fitted between 0.5-100 keV, and Model M was fitted between 0.6-100 
keV since we found strong residuals near ~0.5 keV for this model. The lower energy data (< 2 keV) were modeled using a soft power-law 
(PL) and Gaussian lines to simulate photoionization by the AGN. The top panels of each plot show the unfolded model in EFe units, 
while the bottom panels show the fit residuals in terms of sigma with error bars of size one. Golor scheme: black (NuSTAR PPM A), red 
(NuSTAR FPMB), blue (XMM-Newton PN), green (XMM-Newton MOS), purple (Chandra), cyan (Swift-BAT). 


Table 2 

X-ray spectral fitting results for the AGN in NGC 5643. 


Component 

Parameter 

Units 

Model P 

Model T 

Model M 

Absorber/Reflector 

Nnieq) 

lO^"* cm-2 

- 

lOOt^s 

lo.or^s 


7Vh(1os) 

10^4 cm-2 

- 

1001^5 

5.8t“2 


^inc 

deg 



65.9l^ 


^tor 

deg 

- 

60.0t°;“i 

- 

AGN Continuum 

Dsoft 


3.19±0.13 

q 9Q-I-0.14 
O.ZO_o.i2 

o.04_q 25 


Phard 


1.79±0.05 

1 q-7-1-0.03 
i-y /_0.05 



-f'0.5-2,obs 

10^® erg 

0.4 

0.4 

0.4 


^2-10,obs 

10^° erg 

1.7 

1.7 

1.6 


^30-100,obs 

10^° erg 

15.9 

16.2 

16.3 


i'0.S-2,int 

10^^ erg 

0.4 

0.6 

0.3 


i2-10,int 

10^^ erg 

1.7 

0.8 

0.9 


i30-100,int 

10^^ erg 

1.6 

0.8 

0.6 

r*XMM 

^NuSTAR 



0.94±0.05 

_o.o4 

0.94±0.05 

^Chandra 

^NuSTAR 



0.91±0.10 

u-y4_o.o9 

0 QQ+U J-U 

u.yu_o.o5 

r*BAT 

^NuSTAR 



-I Qq-|-U.J2 

-‘-•^0-0.31 

1.20±0.30 

1 1 c;+0-28 

Xr/d.O.L 



1.10/500 

1.16/499 

1.21/471 


Note. — fixed, ’^unconstrained. Best-fitting model parameter s for Model P (PEXRAV mo del bv IMagdziarz Zdziarskilll995l) , T (TORUS 
model bv IBrightman &: Nandrall2QllJl and M (MYTORUS model bv IMurphv &: YaQOobll2009l) . Details of each model are described in Section 
3. 
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3.1. Model P 

In our fitting of the AGN broadband spectrum, we 
first consider the pexrav r eflection model (‘Model P’; 
iMagdziarz fc Zdziarskil[r995l) . This model has commonly 
been used to model reflection-dominated spectra. How¬ 
ever, it assumes reflection off a slab geometry with an 
infinite column density, which is unlikely to represent 
the true geometry of the AGN torus. It also does not 
incorporate fluoerescence emission lines expected from a 
GTAGN such as the Fe Ka and K/3 lines. Because of 
these limitations, we also fitted the spectra of the AGN 
using more physically motivated reflection/obscuration 
models as described in Section 3.2 and 3.3. However, 
Model P can provide useful comparison with the other 
reflection models and previous studies of this source, as 
well as other CT sources. 

We started by fixing all the line energies and normal¬ 
izations of the soft Gaussian components to the values 
presented in MI3. We also fixed the redshift at z = 
0.0040. Initially, we set the intrinsic line widths to cr = 
50 eV, as we expect them to be unresolved. We then 
allowed the widths for each line to vary in turn, to im¬ 
prove the fit. For those that returned values deviating 
significantly from 50 eV, we fixed them to this new value. 

We added several more Gaussian components to ac¬ 
count for the Fe Ka and Fe K/3 lines, the Ni Ka line and 
also the Compton shoulder associated with the Fe Ka 
line which are not included in the modelH The parame¬ 
ters were also fixed to values obtained in MI3, except for 
the line energy and normalization of the Fe Ka line which 
were allowed to vary throughout. We fixed the inclina¬ 
tion of the reflector to = 65 °. which is equal to that 
modeled bv iFischer et al.l (j20I3[ l for NGC 5643 (with an 
uncertainty of ±5°) based on mapping the kinematics 
of the Ha and [Oiil] narrow line region (NLR) emission 
observed by the Hubble Space Telescope (HST). The re¬ 
flection scaling factor was fixed to i? = — 1 to simulate 
a reflection-dominated spectrum. Model P provides a 
good fit to the data (x^/d.o.f = 551/500), indicating that 
the spectra can be well-fitted by a pure reflection model 
without the need for any direct component. 

The cross-calibration constants of XMM-Newton, 
Chandra and Swift-BAT relative to NuSTAR are con¬ 
sistent with each other with C 1 within the statistical 
errors (see Table 2), indicating that the spectra are in 
good agreement with each other (see also iMadsen et aP 
for the current cross-calibration status between 
different X-ray observatories and instruments with 
respect to NuSTAR). We estimated the intrinsic X-ray 
luminosities from this model by assuming that the 
observed luminosities at 0.5-10 keV and 30-100 keV 
are about 1% and 10% of the intrinsic luminosities at 
these energy bands, respectively (see Table 2). These 
inferred fractions were taken from detailed studies of 
NGC 1068 and have been widely used to estimate the 
intri nsic luminosities of other reflection-dom inated AGN 
(e.g. iMatt et mil997t iBalokovic et al.ll2014ll . 


^ We note that modeling the Compton shoulder using a Gaussian 
component is incorrect since it is not exactly Gaussian, but we 
adopted this approach as an approximate parameterization for this 
particular model. 


3.2. Model T 

We then proceeded to fit the AGN broadband spec¬ 
trum using more physically mot ivated models, start¬ 
ing w ith the TORUS model by iBrightman fc Nandral 
(|2011h (Model T). This model simulates obscuration by 
a medium with a conical section cut from a sphere, with 
a variable biconical polar opening angle (0tor) ranging 
between 26-84°. Small values of 0toi correspond to a ge¬ 
ometrically thick torus, while large values indicate a geo¬ 
metrically thin torus. The line-of-sight A^h(Ios) through 
the torus, which is equal to the equatorial column den¬ 
sity 7VH(eq), extends up to 10^® cm“^, allowing inves¬ 
tigations of the most extreme obscuration, and is de¬ 
fined such that it is independent of the inclination an¬ 
gle (0inc). The model is defined between 0.1-320 keV 
and self-consistently predicts lines commonly found in 
obscured AGN, such as Fe Ka, Fe K^ and also Ka emis¬ 
sion from several other elements (G, O, Ne, Mg, Si, Ar, 
Ga, Gr, Fe and Ni). 

The soft components of the spectrum (< 2 keV) were 
modeled as described earlier in Section 3 and 3.1. Since 
0 inc and 0tor could n ot be constra i ned s imultaneously, 
we fixed 0inc to 65° (iFischer et al.ll20i^ . This model 
also yields a good fit to the data (x^/d.o.f = 579/499) 
with most ht residuals lying at the soft energies (< 2 
keV) and around the iron line emission (~6-7.5 keV; e.g. 
X^/d.o.f = 447/417 when fitting the model by ignoring 
the data around the iron line emission). The fitted intrin¬ 
sic continuum power-law emission from this model has a 
photon index of F = 1.97lgn5 i® absorbed by a col¬ 
umn density of Ah(Ios) > 1.5 x 10^® cm“^; i.e., heavily 
Gompton-thick. The upper limit is unconstrained, with 
values of up to 10^® cm“^ allowed by the model. We 
found that the best-fit 9tor measured by this model is « 
60°. As with Model P, the cross-calibration constants of 
each observation with respect to NuSTAR are consistent 
with 1. We inferred the intrinsic luminosities from this 
model based upon the best-fit parameters obtained, as 
presented in Table 2. 

We also tried to fit the spectrum using different Oi^c 
to see how the key parameters would change. Setting 
0inc to a lower value of 45° produced an unacceptable fit 
(x^/d.o.f = 3989/499). Fixing to a higher value of 
85° to simulate a near edge-on torus inclination, how¬ 
ever, provides a slight improvement to the fit (x^/d.o.f 
= 549/499). The column density obtained is consistent 
with what was inferred before, but is more constrained, 
iVH(los) = 1.8to'6 ^ cm“^. The photon index is 

slightly lower, F = 1.75lg;o3, and the opening angle of the 
torus suggests a surprisingly thin torus, 0tor = 78.5to'74 
degrees. However, we note that iLiii im dMll claim 
that for an edge-on torus inclination, this model overes¬ 
timates the reflection component, and therefore the pa¬ 
rameters obtained may not be reliable. We thus favor 
the best-htting model solution using Oi^c = 65°. 

A closer look at the fitted spectrum between 6-7.5 keV 
(Figure 3) indicates that the model overpredicts the data 
at Ri6.3 keV, which encompass the Gompton shoulder as¬ 
sociated with the Fe Ka line emission, and underpredicts 
the data at Ri6.4-6.8 keV, which include the Fe Ka line 
emission. We investigated whether the observed resid¬ 
uals could be caused by the presence of emission lines 
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Figure 3. A zoom in of the spectral fit for Model T between 
6—7.5 keV showing the asymmetric residuals around the iron line 
emission; the x-axis has been de-redshifted to show the rest-frame 
energy. The color scheme of the data and the model are the same 
as Figure 2. The residuals are dominated by the XMM-Newton 
data at ^6.3-6.8 keV, with the model overestimating the emission 
at 5^6.3 keV, and underpredicting the data at p^ 6.4-6.8 keV. These 
residuals are diminished when applying an ~ H-ll eV offset in the 
XMM-Newton PN data. 

at 6.70 and 6.96 keV corresponding to ionized iron by 
adding Gaussian components at these energies. However, 
we found that the normalizations of the components are 
consistent with zero. We also tried to add a line smooth¬ 
ing component to the model (gsmooth; energy index, 
a = 1) to see if the statistics could be further improved, 
but the fit indicated that this is not required (o-gsmooth ^ 
40 eV). We noticed that the residuals are dominated by 
the XMM-Newton data, and therefore checked whether 
these could be caused by a shift in the energy scale of 
the data. Indeed we found that, applying an ^ -|-11 eV 
offset in the PN data can improve the overall quality of 
the fit (x^/d.o.f = 551/497) and dimi nish the residuals 
around the iron line complex (see also lBauer et al.l[2014l 
for similar finding). However, this energy shift is consis¬ 
tent with the current calibration uncertainty of XMM- 
JVewtor0, and the final results of the fitting are consistent 
with those presented in Table 2. Therefore, we decided 
not to apply this energy shift to the XMM-Newton PN 
data in the final htting of all models. 

3.3. Model M 

We next fitted the AGN spec t rum u sing the MYTORUS 
model by iMurohv fc Yaaoo^ (|2009H (Model M). This 
model simulates obscuration toward an AGN using a 
toroidal absorber geometry (circular cross-section) with 
a fixed opening angle 0tor = 60°. The advantage of this 
model over Model T is that it allows the separation of 
the direct (mytz), scattered (myts) and line emission 
(mytl) components, allowing more freedom to explore 
complex geometries in the modeling. On the other hand, 
the line-of-sight Vh(Ios) is tied to the inclination angle 
and the equatorial VH(eq), which is available only up to 
an absorbing column density of VH(eq) = 10^® cm“^. 
The MYTL component self-consistently includes neutral 
Fe Ka and Fe K/3 fluorescence lines, and the associated 

® The current calibration documentation of XMM- 
Newton EPIC cameras can be downloaded from 
http: //xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 


Gompton shoulders. The model is defined between 0.5- 
500 keV, but we noticed strong residuals in our fit near 
the lower energy threshold, and therefore restricted our 
fits to above 0.6 keV for this model0 

We fitted the AGN spectrum using the simplest ver¬ 
sion of the MYTorus model, which couples all the pa¬ 
rameters of the scattered and fluorescence line compo¬ 
nents to the direct continuum component. We added 
an extra Gaussian component to account for the Ni Ka 
line which is not included in the model (parameters were 
fixed to values obtained in M13). The relative normaliza¬ 
tions of MYTS (As) and mytl (As) with respect to mytz 
(Az) were set to 1. Model M also fits the data well with 
X^/d.o.f = 569/471. The best-fitting global column den¬ 
sity inferred is VH(eq) = 10 . 0 / 1 q 5 x 10 ^^ cm“^, which is 
at the upper limit of the MYTORUS model. Reassuringly, 
the best-fitting inclination angle is 0inc = 65.9~*~? g? de¬ 
grees, consistent with that determined by iFischer et al.l 
(|2013h . However, we note that the MYTORUS model does 
not allow for complete freedom in the fitting of the incli¬ 
nation angles. Due to the assumption made for dtor, the 
model treats a torus with dine < 60° as unobscured, and 
as a result, the best-fit dine measured by the model for 
obscured AGN are usually in the range of ~60-70°. The 
Ah(Ios) measured from this model is well within the CT 
regime, with an unconstrained upper limit; i.e., Yh(Ios) 
= 5 . 8 / 1“2 xlO^^ cm“^. Similar to Model T, the residu¬ 
als of the fit are dominated by the soft energy emission, 
and the emission around the iron line, which can be im¬ 
proved if an energy offset is allowed in the XMM-Newton 
PN data (see Section 3.2). All parameters obtained by 
Model M, including the cross-calibration normalization 
constants and intrinsic luminosities of the AGN, agree 
very well with Model T (see Table 2). 

3.4. The SPHERE Model 

I n addition to the TORUS mo del described in Section 
3.2. lBrightman fc Nandrr] (j2011tl also present a model in 
which the source is fully covered by a spherical geometry 
(^tor = 0°) with variable elemental and iron abundances 
with respect to hydrogen. This model is referred to as 
the SPHERE model. We tried to fit the spectra using 
this model. However, the fit was very poor (y^/d.o.f = 
1169/501) and difficult to constrain. Leaving the metal 
abundances free to vary, which is allowed in this model, 
did improve the fit slightly, but the fit was still poor 
(x^/d.o.f = 848/499), and returned very low values of 
Vh and F, and a high Fe abundance (^ 7). Therefore, 
we will not discuss this model any further. 

4. NGC 5643 X-1 

In addition to the AGN, we also analyzed the NuSTAR 
data for the ULX candidate, NGC 5643 X-1. ULXs are 
off-nuclear point sources with X-ray luminosities exceed¬ 
ing the Eddington limit for the typical ~1OM0 stellar- 
remnant bla ck holes obse rved in Galactic black hole bi¬ 
naries (e.g. [OroiillOOl); i.e.. Lx > 10^® erg s"!, po- 
tentially implying exo tic super-Eddington accretion (e.g. 
iPoutanen et al.ir20Q7[) . While the majority of ULXs haye 

^ Since the fit for Model M was restricted to E > 0.6 keV, a 
Gaussian component that was added to model the emission line 
at E ^ 0.58 keV (M13) in Model P and T was excluded for this 
model. 
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Figure 4. Left: X-ray luminosity variations of NGC 5643 X-1 in the 3-8 keV band observed between 2003—2014 by XMM-Newton, 
Chandra, Swift-XRT and NuSTAR. Observations marked with stars are those that we used to fit the broadband spectrum of the source. 
Right: Best-fit diskpbb model to the combined Swift-XRT-\-XMM-\-NuSTAR spectra of NGC 5643 X-1. The top panel shows the unfolded 
model in EFe units, while the bottom panel shows the fit residuals in terms of sigma with error bars of size one. The data have been rebinned 
to a minimum of 3-sigma significance with a maximum of 50 bins for visual clarity. Color scheme: red (NuSTAR), black (Swift-XRT), 
green (XMM-Newton MOS). 


Table 3 

X-ray spectral fitting results for NGC 5643 X-1. 


Parameter 



Unit 

SwiH-XRT+NuSTAR 

SwiH-XRT+XMM+NuSTAR 

Model = TBABS 

X 

ZPOWERLAW 




A^h 




< 23.2 

< 0.46 

r 




9 

z.oo_Q 33 

1 yn+O.OQ 

4. Qg 

To. 5-24 

r'XMM 

^NuSTAR 



10^° erg 

t.-i-O-o.si 

-I 99 -I-O.O 7 

4.^^-0.03 
-I q-| -|-0.31 
-‘^•'^-‘^-0.23 





Xr/d.O.L 




1.32/7 

1.08/92 

Model = TBABS 

X 

CUTOFFPOWERLAW 




A^h 



10^^ cm-^ 

< 18.5 

< 0.22 

r 




1 70+^-^® 
^•'^-1.64 

1 £ 9 - 1 - 0.13 

4-OZ-O 11 

Ecut 



keV 

10.7l“g 

40.0_g g 

To. 5-24 

r'XMM 

^NuSTAR 



10^*^ erg 

1.08l“;== 

1.18±0.07 

1 1 9 -I-O .36 
-^•-^^- 0.20 





Xr/d.O.L 




1.41/6 

1.03/91 

Model = TBABS 

X 

DISKPBB 




AfH 



10^^ cm-^ 

< 16.8 


Tin 



keV 

4 02^^’*^ 

^•^^-1.47 

^•^^-0.20 

P 




< 0.74 

0.55±0.01 

To. 5-24 

/-•XMM 

^NuSTAR 



10^*^ erg 

1 qi+0-36 

-L-Di_o.20 

1 Qy-l-0.06 

4.U/_o_o8 

1 10-1-0.24 

4-4y_o.23 





Xr/d.O.f. 




1.28/6 

0.99/91 


Note. — ’^unconstrained. 


luminosities of the order o f 10^® erg s“^ (|Walton et al.l 
l2011bl : iSwartz et al.ll2niTl i. a small subset of the pop- 
ulation has been obs e rved to have Ly > lO^o erg s ^ 
(I Walton et al.l l2nnat I.Tonker et al.l 120121 : Button et al.1 
1201211 . Given their luminosities, it has been suggested 
that these more luminous sources might be good can¬ 
didates for hosting intermediate mass black holes (10^ 
< Mq < 10^; e.g. IMiller et al.ll2004 IStrohmaverl 1200911 
accreting at sub-Eddington rates. With a luminosity of 
-bo.5-10 = 2.6 X 10"^° erg s“^ as measured by G04, NGG 
5643 X-1 would therefore be a member of this latter pop¬ 
ulation. 

We extracted the NuSTAR spectra of NGG 5643 X-1 


by defining a smaller circular region of 20"-radius around 
the source using the AGN as a reference point for its 
relative position. For the background region, we used 
an annular segment centered on the AGN (excluding a 
region of '^50"-radius around NGG 5643 X-1 which ac¬ 
counts for the ^70% encircled energy fraction), and with 
the same radial distance (relative to the AGN position) 
and width as the source region. This background region 
was designed to account for contamination by the AGN. 
We then determined the significance of the NuSTAR de¬ 
tection by calculating a no-source probability assuming 
binomial statistics (P > 1% for non-detection; i.e. < 
2.6a, following other NuSTAR studies of fa i nt sources - 
e.g. ILuo et all 120131: ILansburv et al.l [2014 IStern et "ni 
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I2014fl . While the source is significantly detected in both 
of the NuSTAR observations in the 3-8 keV band, it is 
only significantly detected in the second observation in 
the 8-24 keV band (P = 2.21 x 10“® in the combined 
FPMA+B image). 

We next combined the data from the two FPMs for 
each epoch and binned the spectra to a minimum of 5 
counts per bin due to low counts. We then fitted the 
spectra from each epoch using a simple power-law model, 
absorbed by the Galactic column. The spectral fitting 
parameters were calculated using C-stat istics, appropri¬ 
ate for low count statistics (|Cashl [197^ . Although the 
source is not formally detected in the 8-24 keV band in 
the first epoch, we modeled the spectrum including the 
data point in this band to provide a better constraint 
to the fit. We found that the spectra of the source dur¬ 
ing the first and second epoch are consistent with each 
other, F = 2.?;°;® and F = 1.71°;®, / 3-24 = 3.61®;® x 
10“^® erg s“® cm“^ and 3.8^2'q x 10“®® erg s“® cm“^, 
respectively, suggesting that the source has not varied 
significantly between the two observations. The lack of 
a detection at 8-24 keV in the first epoch may be due 
to the location of NGC 5643 X-1 relative to the NuS¬ 
TAR detector gap (Figure 1), which reduced the overall 
sensitivity of the data. 

Motivated by the lack of significant flux variability, we 
therefore combined the NuSTAR spectra of NGG 5643 
X-1 from the two epochs and fitted the total spectrum 
with the Swift-XRT data that were taken simultaneously 
with the first epoch. We reduced the Swift-XRT data as 
detailed in Section 2.4, and extracted the spectrum using 
a 20"-radius circular source region to match the NuSTAR 
extraction region. The background photons were selected 
from a source-free 100"-radius circular region close to the 
source. We grouped the data to a minimum of 40 and 20 
counts per bin for NuSTAR and Swift-XRT, respectively, 
allowing the use of statistics. 

We first fitted the broadband spectrum of NGC 5643 
X-1 over 0.5-24 keV using a simple power-law model ab¬ 
sorbed by a column density intrinsic to the host galaxy 
(tbabs) in addition to the Galactic absorption. We fixed 
the cross-normalization constants between NuSTAR and 
Swift-XRT to 1. The fit is acceptable with y^/d.o.f 
= 9.21/7. The spectrum seems to drop off at E « 11 
keV, which could be an indication of a spectral cut-off 
at around 10 keV, as found for ot h er JVuSTAR-observed 
ULXs fe.g. iWa lton et a l.l 120131: iBachetti et al.l 120131 
IWalton et al.ll20l4lRana et al.ll2015[r ' If shown to be sta¬ 
tistically significant, this would be the first time that such 
a cut-off is observed for NGC 5643 X-1. We therefore 
proceeded to fit the spectrum using an absorbed cut-off 
power-law model. The fit is also acceptable (x^/d.o.f = 
8.46/6), and returned a cut-off energy of Ecut = 10.7/^“ g 
keV. We also fitted the spectra using a multi-color black- 
body accretion disk m odel with a variable te mperature 
disk profile, diskpbb (jMineshige et al.l 11994( 1. which is 
commonly used to model ULX spectra. This model pro¬ 
vides a marginal improvement to the fit over the simple 
and cut-off power-law model (x^/d.o.f = 7.67/6). The 
measured temperature profile, p, i s less than that ex¬ 
pecte d for a thin disk (p < 0.75; iShakura fc Sunva^ 
[1973), consistent with an accretion disk in which advec- 
tion of radiation is important, as expected at very high 


accretion rates where radiation press ure dominates and 
modi fies the structure of the disk (e.g. lAbramowicz et ^ 
Il988f) . This would be consistent with NGC 5643 X-1 ex¬ 
hibiting high-Eddington rate accretion onto a stellar rem¬ 
nant black hole. The luminosity inferred by this model is 
Ao.5-24 = 1 -01-O.2O ^ 10^® erg s“®, assuming that NGC 
5643 X-1 is at the distance of NGC 5643. 

To obtain a better constraint to the broadband spec¬ 
tral fit of NGC 5643 X-1, we fitted the Swift-XRT and 
NuSTAR spectra simultaneously with the 2009 XMM- 
Newton spectrum of the source, in which M13 measured a 
luminosity of L 2-10 = (8.11 ± 0.23) x 10®® erg s“®, con¬ 
sistent with that inferred by our Swift-XRTNuSTAR 
best-fit DISKPBB model L 2-10 = 7.53l^;3g x 10®® ergs“®. 
Data were reduced as described in Section 2.2. Because 
NGC 5643 X-1 fell on the detector gap in the PN camera, 
we only extracted the spectra of the source from the two 
MOS cameras, and co-added the spectra together. We 
extracted the spectra using a 20"-radius circular source 
region, and a 45"-radius circular background region. We 
left the cross-normalization between the XMM-Newton 
spectrum relative to NuSTAR, free to vary. 

We found that both the cut-off power-law and 
DISKPBB models provide better fits to the Swift- 
XRT-\-XMM-\-NuSTAR data than a simple power-law 
model (x®/d.o.f = 93.6/91, 90.2/91 and 99.0/92, respec¬ 
tively). The E-test null hypothesis probability between 
the simple power-law and cut-off power-law model is 
0.023. This low value is strongly suggestive of the need 
for the high energy cut-off, though deeper observation 
will be required to validate this. measured by 

the two models are, within the uncertainties, consistent 
with 1, providing further evidence that the spectra are in 
good agreement with each other. The cut-off power-law 
model measured a cut-off energy of Ecut = 16.3/^9®g® keV, 
and the diskpbb model inferred a temperature disk pro¬ 
file, p = 0.55±0.01 and an inner disk temperature, Tin = 
4.98l^'gg keV. This temperature is consistent with that 
measured for some ULXs when fitting their whole broad- 
band spectrum with just a single dis k component (e.g. 
IWalton et al.l 1201^ iRana et alT 201511 . The best-fitting 
parameters for all models are presented in Table 3, and 
the broadband Swift-XRT-\-XMM-\-NuSTAR spectrum 
of NGC 5643 X-1 fitted by the best-fit model, namely 
DISKPBB, is shown in Figure 4. 

Strong X-ray luminosity variations are often observed 
in ULXs. Motivated by this, we also reduced the 2003 
XMM-Newton data (ObsID 0140950101; UT 2003-02- 
08), as well as the Chandra data, and fitted the 3-8 
keV spectra of NGC 5643 X-1 using a simple power- 
law model absorbed by the Galactic column, to infer its 
3-8 keV luminosities at these epochs. In addition to 
these observations, the galaxy has also been observed by 
Swift-XRT between 2008-2013 on ~30 occasions (mostly 
^1 ks exposure time), with most observations close to 
each other in 2013 targetted on an ongoing supernova in 
the galaxy. For simplicity, for each year in which there 
are multiple Swift-XRT observations of the source, we 
only analyzed the data with the highest nominal expo¬ 
sure time to provide the best estimate of NGC 5643 X-1 
luminosity in that particular year (ObsID 00037275001, 
00037275002, 00037275004 and 00032724009; UT 2008- 
03-16, 2009-09-20, 2010-10-01 and 2013-02-26, respec- 
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lively). Among these observations, NGC 5643 X-1 is 
only significantly detected by Swift-XRT in the 2013 ob¬ 
servation. Therefore for this observation, we extracted 
the data as described earlier, and fitted the 3-8 keV spec¬ 
trum using a simple power-law model absorbed by the 
Galactic column. For the other observations where NGC 
5643 X-1 is not detected, we provided the upper limit lu¬ 
minosities at 90% confidence, estimated using aperture 
photometry assuming F = 2. 

We compare the 3-8 keV luminosities between all the 
observations described above, as well as the NuSTAR 
observations at each epoch in Figure 4 to show the long 
term luminosity variability of NGC 5643 X-1. Note that 
the photon indices were left free to vary in all spectral 
fittings. In general, the plot shows that NGC 5643 X-1 
can vary by a factor of ^2-3 (up to ~5) between epochs. 
This is broadly simi lar to the level of variability observed 
in other ULXs (e.g. iKong et '^11201(11 : iSntton et al.ll2ni3l : 
IWalton et al.ll201^ . These results, combined with the 
parameters inferred by the diskpbb model and evidence 
of an energy cut-off at E ^ 10 keV, as well as the lack of a 
counterpart at other wavelengths (G04), further support 
the ULX classification of NGC 5643 X-1. 

5. DISCUSSION 

In Section 3, we first investigated NuSTAR obser¬ 
vations of the CTAGN candidate in NGC 5643. We 
combined our new data with archival Chandra, XMM- 
Newton and Swift-BAT data, and performed a broad¬ 
band (~0.5-100 keV) spectral analysis of the AGN. On 
the basis of three different reflection/obscuration mod¬ 
els, we found that the AGN is consistent with being CT 
with a column density of A^h(Ios) > 5 x 10^"^ cm“^. The 
absorption-corrected 2-10 keV luminosity obtained from 
the models range between ^ 2-10 = (0.8-1.7) x 10"^^ erg 
s“^, although we note that this could be a factor of a few 
lower/higher due to the fact that the spectra are domi¬ 
nated by reflection with negligible contribution from the 
direct emission, causing large uncertainties in the mea¬ 
surements of the intrinsic emission from the AGN. The 
intrinsic luminosity inferred is at the lower end of the 
luminosity range of t he local bona-fide CT AGN popula¬ 
tion (see Figure 4 in iGandhi et al.l [201411 . Despite the 
presence of a nearby ULX with a 2-10 keV luminosity 
comparable to the observed luminosity of the AGN itself 
(G04), the spectrum of the ULX drops off at U ~ 10 keV, 
and the overall emission above this energy is dominated 
by the AGN. 

This paper provides the first intrinsic X-ray lumi¬ 
nosity measurements for the AGN in NGC 5643 using 
Compton-thick solutions. To provide a basic test of 
the intrinsic luminosity estimate from our X-ray spec¬ 
tral fitting, we can compare the luminosity obtained to 
predictions from multiwavelength indicators. The ob¬ 
scuring torus absorbs disk radiation from the nucleus, 
but then re-emits it in the MIR. We can therefore use 
the intrinsic MIR:X- rav luminosity relationship found by 
IGandhi et al.1 (j2009H based upon high angular resolution 
MIR I2/rm observations to estimate the intrinsic X-ray 
luminosity of the AGN. The advantage of such obser¬ 
vations is that because of their high angular resolution 
(^0.4"), they provide MIR fluxes that are intrinsic to the 
nucleus, minimizing contamination from the host galaxy. 
The 12/im nuclear luminosity for NGC 5643, measured 
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Figure 5. Top: Chandra 6—7.5 keV images of the AGN (left) and 
the simulated PSF (right). The green line shows the direction of 
the semi-major axis of the simulated PSF that was used to pro¬ 
duce the ID radial profiles for both the source and the simulated 
PSF. Bottom: The ID radial profiles of the source (red) and the 
simulated PSF (blue) normalized to their respective total counts. 
The error bars indicate the 1-cr poisson errors for the data. 


from VLT VISIR and Gemini T-R eCS observations, i s 
— ( 1.5T0.4) X 10 " ^^ erg s ^ (|Asmus et al.ll201^ . 
Using the IGandhi et al.l (|2009ll relationship, this corre¬ 
sponds to an X-ray luminosity of T 2-10 = (0.6-2.5) x 
10^^ erg s“^, which agrees well with the 2-10 keV intrin¬ 
sic luminosity inferred from our broadband X-ray spec¬ 
tral modeling^ 

We also compared our intrinsic 2-10 keV luminosity 
with that predicted from the optical [Oiii]A 5007A emis¬ 
sion line luminosity. [Oiii]A 5007A in AGN are mostly 
produced in the NLR by photoionizing radiation from 
the central source. Although this line is not affected 
by the torus obscuration, it can still suffer from ob¬ 
scuration by the host galaxy, which can be corrected 
for using the Balmer decrement. We obtained our 
correc ted [Oiii]A 5007A luminosity from iBassani et al.l 
(|1999ll . T[oiii] = (1.6±0.2) x 10^^ erg s“^. Using the 
L 2 -io:-b[oiii] relationship from iPanessa et all (I2006D . we 
infer an X-ray luminosity of T 2-10 = (2.1-33.9) x 10^^ 
erg s“^The lower end of this luminosity range is consis¬ 
tent with the luminosity measured from our X-ray spec- 

The given luminosity range accounts for the mean scatter of 
the correlation. 
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Figure 6. Optical image of NGC 5643 retrieved from the ESO 
Digital Sky Survey image server. Chandra position of the AGN 
(NGC 5643) and the ULX candidate (NGC 5643 X-1) are both 
labeled and marked with an “x”. North is up and East is to the 
left in the image. 

tral fitting. This is supported by the MIR [Oiv] A25.89/rm 
emission line, Wl = (2.7±0.1) X lO^o erg s'^ 
(|Goulding fc Alexander] [200^ . Using L2-io:7y[oiv] I'ela- 
tionship defined by iGonlding et al.l (j201flrF1 . we deter¬ 
mine an X-ray luminosity of T 2-10 = (0.9-5.5) x 10"^^ 
erg s“^.® This luminosity agrees very well with that mea¬ 
sured by our best-fit models. 

With our updated intrinsic X-ray luminosity, we can 
estimate the Eddington ratio of the AGN. Using Mbh 
estimated from the M bh-ct* relation, Mbh = Mq 
(|G oulding et al.l [2010t see also footnote 3), combined 
with the range of intrinsic X-ray power that we measured 
from the best-fit models, we infer that the AGN is ac¬ 
creting matter at about 5-10% of the Eddington rate.^^ 
However, we note that this value is subject to large errors 
(factor of a few) due to th e highly uncertain quantiti es 
involved in the calculation ([Brandt fc Alexandedl201^ . 

In the past few years, much evidence for the presence of 
distant cold r eflecting regions have been found in nearby 
GTAGN le.g. iMarinucci et al.ll201^ lArevalo et al.l[20l4 

iBauer et al.ll2014h . To accurately characterize the broad¬ 

band spectra of these AGN, multiple reflector compo¬ 
nents are required to model each of these reflecting re¬ 
gions. We investigated whether this could be the case 
for NGG 5643 by looking at the Chandra image of the 
AGN in the 6-7.5 keV band (^30 counts). This is shown 
in the top left panel of Figure 5. At first glance, the 
image appears to show that the iron emission region is 
extended well beyond the parsec-scale torus, suggesting 
the presence of distant cold reflecting regions. However, 
this appearance could also be due to an elongated PSF 
as a result of the position of the target on the detector, 
as well as the energy range used to produce the source 
image. We therefore used the MARX simulatoiQ (v4.4) 
to produce the expected PSF with a much higher count 
than the source counts (~ 100 000 counts) at the po- 

We assume a typical AGN bolometric correction L^^ilL 2—10 

20 to calculate the bolomet ric luminosity and hence the mass 
accre tion rate of the AGN (e.g. lElvis et al.lll993 : IVasudevan et al.l 

[20ICT) . 

The MARX simulation documentation is available at 
http: //space.mit.edu/ASG/MARX/ 


sition of the AGN, and at the mean energy range that 
we are interested in {E = 6.75 keV). This is shown in 
the top right panel of Figure 5. We then plotted ID ra¬ 
dial profiles of the target and the simulated PSF along 
the semi-major axis of the simulated PSF determined us¬ 
ing the WAVDETECT tool within CIAO, and performed a 
Kolmogorov-Smirnov (KS) test between the two distri¬ 
butions (bottom panel of Figure 5). We found that the 
KS test probability is Pks = 0.983, indicating that the 
iron emission region of the target is consistent with the 
simulated PSF, and therefo re is not extended (see simi¬ 
lar finding for NGG 3393 in IKoss et al.ll20I5ll . A deeper 
(scheduled) Chandra observation could help to confirm 
this. 

The GTAGN in NGG 5643 is located in a face-on 
grand design spiral galaxy with an undisturbed disk 
(see Figure 6). The emission lines from the AGN are 
not obscured/extinguished by the host galaxy material 
([Goulding fc Alexanded [200^ . indicating that any phe¬ 
nomenon that might contribute to the GT obscuration 
in this source is therefore likely to occur very close to 
the nucleus. Indeed, there have been indications that 
the emission at the west side of the nucleus is obscured 
by a warped disk or the presence of a dust lane (e.g. 
Simpson et al.l II9971 iFischer et al.l I20I3L iDavies et al.l 
20I4f) . The torus itself could be the inner structure of 
this disk. Because of the more direct and “cleaner” view 
of the central engine in NGG 5643 as compared to some 
other very nearby bona-fide GTAGN such as NGG 1068, 
NGG 4945 and Gircinus, this makes NGC 5643 a good 
candidate for more detailed studies of the CT obscurer 
in the nuclear region. 

6. SUMMARY 

We observed the AGN and ULX candidate in NGC 
5643 using NuSTAR on two occasions conducted at 
about a month separation. A summary of our main re¬ 
sults is as follows: 

1. Using the combined data from NuSTAR, Chandra, 
XMM-Newton and Swift-BAT, we have extended 
the broadband spectral analysis of the GTAGN 
candidate in NGC 5643 to high energies (~0.5-100 
keV). Using physically motivated toroidal obscura¬ 
tion models, we showed that the source is indeed 

CT with a column density of Vh(Ios) > 5 x 10^^ 

_2 
cm . 

2 . The absorption-corrected 2-10 keV luminosity 
measured by these models is T 2 -io,int = (0.8-1.7) 
X 10^^ erg s“^, consistent with that predicted 
from multiwavelength intrinsic luminosity indica¬ 
tors. The luminosity inferred is at the lower end of 
the luminosity range of the local GTAGN popula¬ 
tion. 

3. The NuSTAR spectra of the off-nuclear source, 
NGC 5643 X-1, shows evidence for a spectral cut¬ 
off at A ^ 10 keV, similar to that observed in other 
ULXs observed by NuSTAR. Combining this in¬ 
formation with the evidence for X-ray luminosity 
variations observed between different observations, 
along with the absence of unambiguous counter¬ 
parts at other wavelengths, we concluded that the 







































NuSTAR Observations of NGC 5643. 


13 


source is consistent with being a ULX. Future si¬ 
multaneous low and high energy X-ray observations 
of this held are needed in order to provide higher 
quality data to conhrm the spectral cut-off that we 
observed, and to better characterize the broadband 
spectrum of the source. 
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